Quantum Andreev Oscillations in normal-superconducting-normal nanostructures 
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We show that the voltage drop of specially prepared normal-superconducting-normal nanostruc- 
tures show quantum Andreev oscillations as a function of magnetic field or input current. These 
oscillations are due to the interference of the electron wave function between the normal parts of the 
structure that act as reflective interfaces, i.e. our devices behave as a Fabry-Perot interferometer 
for conduction electrons. The observed oscillations and field periods are well explained by theory. 
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The possibilities of exploiting quantum mechanical ef- 
fects - with all the interferences and other phenomena oc- 
curring in real nano-devices - find new expectations that 
may lead to the fabrication of small devices with appli- 
cations in new fields of technology as ballistic electronics 
and spintronics and flux devices combining normal and 
superconducting materials. Earlier work in semiconduc- 
tors and STM experiments demonstrate the existence of 
quantum oscillations [l], 0] ■ Recently, spin-polarized res- 
onant tunneling in magnetic tunneling junctions showed 
large changes in the magnetoresistance due to the inter- 
ference of the carrier wave function Q . In this work we 
seek after quantum mechanical interference effects in the 
magnetoresistance using normal-superconducting struc- 
tures without tunneling. Because one of the characteris- 
tics needed is ballistic transport our experiments have to 
be done at low temperatures and the systems should be 
designed to have Fermi wavelength of the order or larger 
than their relevant size. 

Assume a strip with a lateral structure M1/M2/M1, 
where Ml and M2 are two different materials with differ- 
ent Fermi energies Ef and lengths L\ and L 2 and where 
the electrical current passes through them. Because of 
the different Ep's between Ml and M2, the one particle 
potential can be described by a barrier t/jj of length L 2 in 
M2 that acts as a potential well where the wave function 
behaves coherently having multireflections, i.e. a kind 
of Fabry-Perot interferometer for electrons, a problem 
studied recently for the case of fluctuations in the mag- 
netoresistance of graphene Q. Inset in Fig. 1(a) shows 
the one-dimensional geometry of the system with the bar- 
rier U 2 {x) depending on the potential drop between the 
two ends of the trilaycr. The transmittivity along such 
structure is 0, [B| 



T(U,a) 



4AE(a)E(a) 



AE{a)AE{a) + C7| sin 2 (27rL JAE{a)^) 



(1) 

The parameters in {T]) are: AE(a) = E{a) — Uz; E{a) — 
{Ef + U) cos 2 (a) is the energy of the incoming particles 



forming an angle a to the interface, m the electronic 
mass and H — h/2n. Applying a potential to the trilayer 
the potential in M2 is U 2 {x) = U - {U/L 2 ){x - L x ) 
{Li ^ x < L\ + Li). The solution to this problem can 
be obtained as combination of Airy functions [J 0] . A 
good approximation is to take the trapezoidal rule, i.e. 
U 2 {x) w U - 0.5U and calculate T{U) following (JTJ as 
shown in Fig. 1(a). Interference effects require ballistic 
transport in the potential region M2 with conveniently 
flat interfaces to avoid multiple reflections. We will use 
M2 as a superconductor (S) or as normal metal (N) and 
vice versa for Ml, of appropriate lengths. 

When M2 is a superconductor the current will be con- 
trolled by Cooper pairs Andreev currents if U < A (A is 
the superconducting gap) 0. However at higher U the 
current is controlled by quasiparticles and for U > 3A 
we have practically conduction between normal materi- 
als @, 0j S| - The solution for the current / is in this 
case: 
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where the term inside the integral is the Andreev conduc- 
tance (in units of quantum of conductance) g^s{U,T) 
0, [H, which depends on T in the normal state; a — 
(U/A) - [(£VA) 2 - I] 5 ; A is a constant that depends 
on the junction geometry and on the integration average 
on a. The solution for / is depicted in Fig. 1(b). Note 
that: (a) g^is oscillates as a function of the energy drop U 
because of the resonances in T(U) and it tends to T(U) 
when the product U/A > 1, see Fig. 1(a). (b) The mea- 
sured voltage drop will oscillate as a function of I and, at 
constant /, as a function of the magnetic field B through 
U{B)/e oc V{B) = IR{B), where R{B) is the, in general, 
non-oscillatory magnetoresistance of the sample. 

The increment of the applied voltage drop AV 
needed to obtain one oscillation is given by AV^ ~ 
it 2 {ti 2 / m) {1 / L 2 ) . It can be seen that for m of order of 
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FIG. 2: (a) Sketch of the superconducting long strip with 
voltages electrodes normal to it. The distance between them 
was 5 /im. The width and thickness were 200 nm and 80 nm. 
The SEM picture at the upper right shows the region at and 
around one of the slits, (b) SEM picture of the WC ring 
labeled WR2 with the contacts pads deposited after the ring, 
(c) A blow out of the same ring as in (b). The cartoon below 
shows the equivalent circuit of the ring with Ml and M2 the 
normal and superconducting paths, (d) SEM picture of ring 
WRo produced in a single-step procedure. The white bars in 
(c) and (d) indicate 300 nm distance. 



FIG. 1: (a) Inset shows the energy diagram for a M1/M2/M1 
structure with a barrier U2(x) such that ^2(^1) = Uq (the 
energy drop at Ml has been neglected for clarity). Main panel 
shows the transmittivity T from JT]) for M2 in the normal 
state (dashed line) as well as the Andreev conductance §ns 
(continuous line) vs. normalized energy drop U/A at M2. 
(b) Carrier current I vs. normalized energy drop after 
The two dotted lines are calculated for constant transmittivity 
T — 0.5, 0.7 for a NS bilayer that shows no oscillations. The 
middle curve that shows the oscillatory behavior is obtained 
with T from (JTJ), see panel (a). Note that the strength of the 
oscillations is larger for U < A. 



the free mass the needed voltage change is AV ~ fiV to 
mV for L2 in the micro and nanometer range, respec- 
tively. For example, if the superconducting gap is less 
that 1 meV with L2 ~ 10 nm we need to have voltage 
drops of several millivolts, i.e. the situation is practically 
a trilayer with normal conducting materials. However, 
for L 2 ~ 1 (mi, AV ~ 1 fxV. 

The NSN and SNS structures were produced combin- 
ing electron beam lithography with the deposition fea- 
tures of a FEI NanoLab XT 200 microscope. The tung- 
sten precursor in the gas injection system enables us 
to fabricate homogeneous, amorphous superconducting 
tungsten-carbide (WC) micro- and nano-structures by 
Ga + -ion beam induced deposition (IBID) with critical 
temperature T c = 4 — 5 K Q. The deposition param- 



eters were 30 kV accelerating voltage, 10 pA Ga + -ion 
current and 50% overlap. All samples were produced on 
5 x 5 x 0.53 mm 3 silicon substrates with 150 nm insulating 
SiN layer. The lithographically prepared multi-electrode 
structure was connected to a multi-electrode chip and 
this was fixed in a magneto-cryostat system. The resis- 
tance measurements were done using four-wires method 
with AC Linear Research LR700 bridge with a multi- 
plexer. DC Current-voltage (I — V) characteristics as 
well as field B dependent V(B, I) curves were obtained 
using a Keithley 2182 nanovoltmeter with a 6221 current 
source. 

We will discuss first the results for SNS-type nanos- 
tructures. A ~ 60 ^m-long WC strip with eight volt- 
age electrodes at a distance of 5 /xm was fabricated in a 
single-step IBID procedure enabling us to measure simul- 
taneously different portions of the strip (channels 1 to 8, 
see Fig. [UJa)). To produce narrow, normal-conducting 
paths (M2) within the superconducting strip with high- 
quality SN interfaces, minimizing preparation time and 
costs, we irradiated it locally at the middle of different 
voltage channels with the same Ga + -ion beam (30 kV, 
1 pA, zero nominal thickness). In general, ion irradia- 
tion strongly degrades the superconducting properties of 
a material, see e.g. Ref. [13. The length of the slits was 
L2 — 30 nm, see Fig. [UJa). As reference some channels 
were left without slits. The sample W-SD-1 discussed 
here showed a T c ~ 4.3 K defined at half-resistance at 
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FIG. 3: (a) Voltage (V) vs. input current (7) normalized by 
the critical current / c (0) = 19.1 (iA at zero field at 3.6 K, 
measured for two different channels, each one with one slit: 
Ch. 2 at T< B ^ 1 T and Ch. 4 at 0.1 T applied field. Two 
curves were obtained for Ch. 8 (without irradiation, I c (0) = 
16 fiA) at T and 8 T of the same strip sample, (b) Voltage 
vs. applied field for Chs. 2 and 8 at 3.5 K and for an input 
current of 12 /xA. 



B = T. At input currents I < 1 fiA and B = 0T the 
slits short circuit the strip at T < 4 K due to the large 
proximity effect as observed, for example, in Al thin films 
At higher currents or fields the slits show normal 
conducting behavior. 

Figure [3Ja) shows the voltage vs. current measured 
for Chs. 2, 4 and 8 at 3.6 K at different magnetic fields 
applied normal to the input current. At zero field and at 
/ < f 6 /iA the voltage drop in the homogeneous, unirra- 
diated Ch. 8 is below the 1 nV resolution. At the critical 
current I c ~ 16 fiA the voltage jumps to a value close 
to that expected for that conducting path in the nor- 
mal state, which has a resistance of ~ 500 ft. Within 
experimental resolution this curve does remain without 
changes for an applied field of 0.01 T. At B > O.Of T 
the critical current decreases and flux-flow behavior is 
observed. As example, Fig.OJa) shows the data obtained 
at B = 1 T. No hysteretic behavior was observed, which 
rules out heating effects. Channels 2 and 4 show the 



voltage across two different SNS structures, each with a 
single slit. At B = T and O.Of T the voltage drop below 
I c (0) is several orders of magnitude larger than that ob- 
tained for Ch. 8 (as example we show the signal of Ch. 2 
in Fig. da)). At fields below f T and at T = 3.6 K Chs. 
2 and 4 show clear oscillatory behavior similar to that 
shown in Fig.QJb), in contrast to that of Ch. 8. At higher 
fields (B > 1 T) the flux-flow contribution prevails and 
the curves for all channels resemble qualitatively (Chs. 2 
and 8 response is shown in Fig.[3]Ja)). The quantum me- 
chanical resonances can be also observed in V{B) at con- 
stant current, see Fig. O^b), because the magnetic field 
changes the voltage drop at the normal conducting path 
due to its intrinsic magnetoresistance changing therefore 
T(U). The oscillation amplitude is of the order of mV, 
i.e. larger than the superconducting gap A < 0.4 meV, 
in agreement with the model if L 2 — 30 nm. Note that 
the free electron gas model used here provides the gen- 
eral physics of the phenomenon. More detailed analysis 
of the band structure of Mf and M2 can explain the del- 
icate structure of the experimental curves. 

In what follows we discus NSN-type nanostructures 
with long superconducting M2 paths. For reasons that 
will become clear below and as an example of the de- 
position possibilities we prepared NSN ring structures, 
see Fig. EJb-d). To do the normal Mf paths we used 
a two-step IBID procedure, which consists in re-deposit 
WC material on the top of the already deposited su- 
perconducting ring path as the electrical contacts of the 
rings WRi and WR2. fn this case the material Mf, be- 
tween the superconducting electrodes and ring, shows 
a lower T c as the resistance transition indicates, see 
Fig. [H a). The T— and B— dependent measurements re- 
veal that we have normal (semi) conducting interfaces be- 
tween the electrodes and the ring structure at T > 4.7 K 
at B = T. To check that the observed two-step tran- 
sition as well as the oscillations described below are not 
related to the ring structure itself, a one-piece (e.g. one- 
step procedure) ring WRo was prepared in which the con- 
tacts and ring are produced in a single-step IBID using 
especial bitmap pattern generator features, see Fig.[2^d). 

Because in our NSN ring structures L 2 ~ 3 fim, we ex- 
pect oscillations in the submicrovolt range, fn agreement 
with the model these oscillations are nicely observed be- 
tween ~ 4.9 K and 4.7 K and around 0.2 T, see Fig.QJb). 
The field was always applied normal to the ring area. No 
oscillations were observed at lower temperatures, below 
the critical line of the Mf material, above the critical line 
of the ring as well as for the homogeneous ring WRo in 
the whole T, jB-range. Note that the model [4[ applies 
only at T = OK. However, the equivalent maximum en- 
ergy of the oscillations in the ring structure is ~ 50 neV 
much smaller than the thermal energy of 0.4 meV. This 
indicates that the conducting path Ml is a narrow band 
semiconductor because only in this case the temperature 
excitation of electrons, intraband or to other conduction 
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FIG. 4: (a) Resistance vs. temperature at different applied 
fields of ring WRi measured with a current of 1 fj,A, which is 
30 times smaller than I c at similar B,T. The critical field of 
the superconducting transition of the ring itself (upper tran- 
sition) follows Bc2(T) ~ 8.1[T](1 - (T/Tc) 2 ) with T c ~ 5.2 K. 
(b) Resistance of the ring WRi vs. field at different T. The 
inset blows out part of the 4.8 K curve measured increas- 
ing and decreasing field. The dimensions of the ring WRi 
were (2.57, 1.6, 0.15) /im for the outer- and inner-diameter 
and thickness, respectively. 



band, is small. 

We expect that the field period of the oscillations 
should be very sensitive to the properties of the normal 
Ml part of the NSN structure and therefore it should 
change from sample to sample. This is indeed the case. 
The ring sample WR2 shows a field frequency of 0.1 T, 
about ten times larger than for sample WRi, implying 
that the Ml material should have a larger Ep according 
to the model. 

The highly reversible behavior of the oscillations at all 
temperatures and fields (an example at 4.8 K is shown 
in the inset of Fig. BJb)) does not support an interpre- 
tation of the oscillations based on Josephson-like junc- 



tion at the contact positions. Magnetoresistance oscilla- 
tions were observed in superconducting mesoscopic rings 
at very low fields and explained by the oscillatory behav- 
ior of the critical current density I C (B) at T ~ T C (B) 
Oscillations in the conductance were also observed 
in normal metal rings with two tunnel junctions due to 
the magneto-electric Aharonov-Bohm effect [l3| ■ In those 
works the field period is related to the field needed to add 
one (or two) flux quantum in the ring area, which for the 
geometry of our rings would give ~ 0.5 mT - 2.5 mT, 
far below the experimental observed value of 12 mT - 
100 mT for the rings WRi and WR2, respectively. 

In summary, in this work we showed the existence of 
quantum Andreev oscillations in the magnetoresistance 
of NSN and SNS nanostructures. These oscillations are 
observed also as a function of magnetic field due to the 
intrinsic, non-oscillatory magnetoresistance effect of the 
structure. The oscillations and their field period are well 
explained by theory taking into account ballistic trans- 
port and carrier wave function interference within the 
structure. The effects here reported might be used to 
study granular, superconducting materials, when neither 
Meissner effect nor percolation can be observed. Simi- 
lar oscillations as in our NSN structures were rece ntly 
observed in thin mesoscopic structures of graphite [14| . 
This kind of quantum mechanical resonances in the volt- 
age should appear in every nanostructure composed by 
different materials, if the electronic conduction is ballistic 
and the roughness of the interfaces is small enough. 
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